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Annihilation of extremely energetic cosmic neutrinos on the relic-neutrino background can give rise 
to absorption lines at energies corresponding to formation of the electroweak gauge boson Z°. The 
positions of the absorption dips are set by the masses of the relic neutrinos. Suitably intense sources 
of extremely energ etic (10 21 - 10 25 -eV) cosmic neutrinos might therefore enable the determination of 
the absolute neutrino masses and the flavor composition of the mass eigenstates. Several factors — 
other than neutrino mass and composition — distort the absorption lines, however. We analyze 
the influence of the time-evolution of the relic-neutrino density and the consequences of neutrino 
decay. We consider the sensitivity of the lineshape to the age and character of extremely energetic 
neutrino sources, and to the thermal history of the Universe, reflected in the expansion rate. We 
take into account Fermi motion arising from the thermal distribution of the relic-neutrino gas. We 
also note the implications of Dirac vs. Majorana relics, and briefly consider unconventional neutrino 
histories. We ask what kinds of external information would enhance the potential of cosmic-neutrino 
absorption spectroscopy, and estimate the sensitivity required to make the technique a reality. 
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I. INTRODUCTION 



A. A Brief History of Relic Neutrinos 

According to the standard cosmology, neutrinos should 
be the most abundant particles in the Universe, after the 
photons of the cosmic microwave background, provided 
that they are stable over cosmological times. 1 Because 
they interact only weakly, neutrinos decoupled when the 
age of the Universe was w 0.1 s and the temperature of 
the Universe was a few MeV. Accordingly, relic neutri- 
nos have been present — as witnesses or participants — for 
landmark events in the history of the Universe: the era 
of big-bang nucleosynthesis, a few minutes into the life 
of the Universe; the decoupling era around 379 000 y [g] , 
when the cosmic microwave background was imprinted 
on the surface of last scattering; and the era of large- 
scale structure formation, when the Universe was only a 
few percent of its current age. 2 

Some of the earliest cosmological bounds on neutrino 



* Gabricla.Barcnboim@uv.es 
^ omena@fnal.gov 
' quigg@fnal.gov 

1 For compact summaries of the canonical thermal history of 
the Universe, see Weinberg [J, §15.6, the review article by 
Steigman 0, §19—23 of the Review of Particle Physics 0, or 
a standard textbook 0000 

2 For a recent quantitative assessment of evidence that neutrinos 
were present at these times, see Ref. 



masses followed from the requirement that massive relic 
neutrinos, present today in the expected numbers, do not 
saturate the critical density of the Universe [iJJ, UJJ • Re- 
fined analyses, incorporating constraints from a suite of 
cosmological measurements, sharpen the bounds on the 
sum of light-neutrino masses 0] . The discovery of neu- 
trino oscillations [l3llT3.IT5j implies that neutrinos have 
mass, but we cannot reliably compute the contribution 
of relic neutrinos to the dark matter of the Universe until 
we establish the absolute scale of neutrino masses. Cur- 
rent estimates for the neutrino fraction of the Universe's 
mass-energy density lie in the range 0.1% < fl u < 1.5%. 

The neutrino gas that we believe permeates the present 
Universe has never been detected directly. Imagina- 
tive schemes have been proposed to record the elas- 
tic scattering of the 1.95-K relic neutrinos, but all ap- 
pear to require significant further technological devel- 
opment before they can approach the needed sensitiv- 
ity [l6Lll7T| . In this paper, we elaborate a complementary 
approach: detecting relic neutrinos by observing the res- 
onant annihilation of extremely-high-cncrgy cosmic neu- 
trinos on the background neutrinos through the reaction 

vv -> Z° mUll^liJIHIa!]. By observin S z " bursts 
or absorption lines, one may hope to determine the ab- 
solute neutrino masses and the flavor composition of the 
neutrino mass eigenstates. 3 



3 See Ref. 1241 for a general review of other aspects of particle 
physics at neutrino observatories. 
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As a Gedankenexperiment, the prospect of cosmic- 
neutrino absorption spectroscopy has great clarity and 
appeal. Reality is more complicated, and it is our 
purpose — building on earlier work — to analyze all the im- 
portant effects that will influence the execution and in- 
terpretation of neutrino-absorption experiments. We are 
encouraged in this effort by the imminent construction 
and operation of neutrino observatories and by imagi- 
native efforts to develop new techniques to detect super- 
high-energy neutrinos. A novel aspect of the analysis pre- 
sented here is our attention to the thermal motion of the 
relics. We also raise the possibility that cosmic-neutrino 
absorption spectroscopy might open a new vista on the 
thermal history of the universe, as well as extending or 
validating our understanding of neutrino properties. 

In the body of this introductory section, we develop the 
pieces that enter the analysis of neutrino absorption spec- 
tra: our expectations for the relic neutrino background 
now and in the past, details of the annihilation cross sec- 
tion, and possible sources of extremely energetic cosmic 
neutrinos. We also survey experiments that aim to de- 
tect ultrahigh-energy neutrinos. In iflTl we describe the 
idealized situation of a super-high-energy neutrino beam 
incident on a (very long) uniform column of relic neutri- 
nos at today's density, but with negligible temperature. 
We describe the information that could be extracted from 
absorption dips, assuming perfect energy resolution and 
flavor tagging. 

The extremely long interaction length for neutrinos 
traversing the relic background means that we must inte- 
grate over cosmic time, or redshift, and this takes up illlll 
There we discuss the mechanisms that distort absorption 
lines and how the distortions compromise the dream of 
determining the absolute neutrino masses. We also re- 
mark on the sensitivity of the line-shape to the thermal 
history of the Universe. 

We include Fermi motion due to the relic-neutrino 
temperature — which evolves with redshift — in i flVI The 
mean relic-neutrino momentum at the present epoch acts 
as a rough lower bound on the effective target mass. Sec- 
tion is devoted to the implications of unconventional 
neutrino histories, including neutrino decay and the con- 
sequences of a lepton asymmetry in the early Universe. 
We summarize what we have learned, and assess the 
prospects for experimental realization of these ideas in 
WII Looking forward to the experiments, we consider 
how external information could enhance the potential 
of cosmic-neutrino absorption spectroscopy, and we es- 
timate the sensitivity required to make the technique a 
reality. 
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where p is the relic momentum and T is the temperature 
of the photon ensemble. The number density of photons 
throughout the Universe is 



rCO = 



1 



d 3 p 



1 



exp (p/T) - 1 



2C(3) T 3 



(2) 



where £(3) « 1.20205 is Riemann's zeta function. In 
the present Universe, with a photon temperature To = 
(2.725 ± 0.002) K @, the photon density is 



n 7 o = n-yiTo) « 410 cm 



(3) 



The present photon density provides a reference for 
other big-bang relics. The essential observation is that 
neutrinos decoupled when the cosmic soup cooled to 
around 1 MeV, so did not share in the energy released 
when electrons and positrons annihilated at T « m e , 
the electron mass. Applying entropy conservation and 
counting interacting degrees of freedom, it follows that 
the ratio of neutrino and photon temperatures (below 
m e ) is 

T u /T=(±) 1/ \ (4) 
so that the present neutrino temperature is 

4 \V3, 



T„o 



T = 1.945 K ~» 1.697 x 10~ 4 eV . (5) 



The momentum distribution of relic neutrinos follows 
the Fermi-Dirac distribution (with zero chemical poten- 
tial), 
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The number distribution of relic neutrinos is therefore 
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In the present Universe, the number density of each (ac- 
tive) neutrino species is 5 



n Ui o = n„c = n Vi (T v o) « 56 cm 



-3 



(8) 



B. Character of the Relic Neutrino Background 

The cosmic microwave background is characterized by 
a Bose-Einstein blackbody distribution of photons (per 



4 We adopt units such that h = 1 = c, and we will mea- 
sure temperature in kclvins or electron volts, as appropriate to 
the situation. The conversion factor is Boltzmann's constant, 
k = 8.617343 x 1CT 5 eV K" 1 . 

5 The unconventional neutrino histories described in Jyjcan alter 
this expectation. 
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FIG. 1: Neutrino temperature T v (right-hand scale) and mean 
momentum p v of relic neutrinos (left-hand scale) as a function 
of the redshift z. 



and the mean momentum of relic neutrinos today is 
7 C(4) 



2 C(3) 



■T v0 w 3.151^0 « 5.314 x 10" 4 eV , (9) 



where we have used £(4) = 7r 4 /90 = 1.08232. In the same 
way, the mean-squared neutrino momentum is given by 

(Plo) = 15 St ■ T% « 12.94 T„ 2 , (10) 



so that 



C(3) 



(p 2 }2 w 3.597 T„ w 6.044 x 10" 3 eV . (11) 



Neutrinos decoupled very early in the history of the 
Universe, at redshift z — O(10 10 ). The temperature of 
a massless decoupled species scales as T ~ (1 + z), as 
shown in Figure ^ f° r redshifts since the formation of 
the first stars By Eq. Q, this means that the relic 
neutrino number densities will be redshifted as n Vi {z) — 
n V Q (1+z) 3 . The dependence of the relic neutrino number 
on the redshift is depicted in Figure [5] 

The effective relic neutrino density that an extremely- 
high-energy neutrino would encounter while traversing 
the expanding Universe is the neutrino density per unit 
redshift. The propagation distance r is related to the 
redshift z through 



dr = dz/(l + z)H(z) , 



(12) 



where H(z) is the Hubble parameter. In a flat Universe 
with negligible radiation component, we may write 



h 2 {z) = hI [n m {\ + zf + n A ] , 



FIG. 2: Relic neutrino number density versus redshift. 



where ft m is the matter density, fl\ is the cosmological 
constant, and Hq = h ■ 100 km/s Mpc" 1 is the present 
value of the Hubble constant. The neutrino density per 
unit redshift, sometimes called the column density, is 
then 



n Vi (z) 



n vQ (1 + zf dr 
n„ (1 + z) 3 dz 

+ W) 



(14) 



n v0 (1 + z) 2 dz 



In Figure |31 we show this column density as a function of 
the redshift for the ACDM model, 6 that is, a flat universe 
with a cosmological constant and cold dark matter. 

The appearance of the Hubble parameter in Eq. (|14|) 
raises the possibility that careful observation of neutrino 
absorption lines could reveal something of the thermal 
history of the Universe. We shall see in illlll examples of 
how the lineshapes differ in different cosmologies. 



C. The Neutrino Mass Spectrum 

Through the past few years, experiments have adduced 
robust evidence for flavor change involving solar, atmo- 
spheric, and reactor neutrinos, as well as neutrinos pro- 
duced by accelerator beams |2g . Putting aside exotic in- 
terpretations, these results establish that neutrinos have 
nonvanishing masses and that neutrino species mix. The 
most economical description of the new phenomena is 
given in terms of the 3x3 neutrino mixing matrix that 



(13) 



6 For three complementary views of today's concordance cosmol- 
ogy, see Refs. lUHH^. 
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FIG. 3: Column density ifTIjl versus redshift for the ACDM 
model, with Q. m = 0.27 ± 0.04, fi A = 0.73 ± 0.04, and h = 



relates flavor eigentstates to mass eigenstates, 7 the ana- 
logue of the (Cabibbo-Kobayashi-Maskawa) quark mix- 
ing matrix. The standard neutrino-mixing phenomenol- 
ogy entails six parameters: three real mixing angles (#12, 
$23, O13), one Dirac CP phase (5), and two independent 
mass-squared differences 8 (Am^ and Am 2 ^). 

To connect the solar, atmospheric, reactor, and ac- 
celerator observations with neutrino parameters, we fol- 
low convention |3l| in identifying the mass splittings and 
mixing angles that drive the solar and atmospheric tran- 
sitions as (Am^ 2 ,6'i2) and (lAm^l, 823). The sign of the 
mass splitting between the atmospheric state V3 and the 
solar doublet {1*1,1/2) is not yet known. Both the nor- 
mal hierarchy (777.3 > 7772 > mi) and the inverted hierar- 
chy (777,3 < mi i5 "72) are illustrated in Figure where 
the colored bars represent the flavor content of the mass 
eigenstates. 

The best fit point for the combined analysis of solar 
neutrino data together with KamLAND reactor data [3i| 
is at Am 2 12 = 8.2+° j x 1CT 5 eV 2 and tan 2 6 l2 = 
0.40+Q Qy. In the atmospheric neutrino sector, the most 
recent analysis of K2K accelerator data and atmospheric 



neutrino data constrains 1.9 x 10 3 eV 2 < |Am| 3 | < 
3.0 x 10~ 3 eV 2 , sin 2 26> 23 > 0.90 [H|35|. 

Only an upper bound exists on the mixing angle 
#13 (which connects the solar and atmospheric neutrino 
realms), and the CP-violating phase 8 is unobservable in 
current neutrino oscillation experiments, so we allow in 
Figure El for the variation -1 < cos 6 < +1. The CHOOZ 
reactor experiment bounds sin 2 2^13 < 0.1 (at 90% CL) 
for a value of the atmospheric mass gap close to the cur- 
rent central value [36|. 

Several oscillation experiments that exploit neutrino 
beams from nuclear reactors and accelerators are taking 
data, and similar experiments will take data over the next 
few years. In particular, future reactor neutrino experi- 
ments could set the value of #i 3 , as explored in detail in 
Ref. [I]} • For a recent study of the measurement of lep- 
tonic CP violation and the pattern of the neutrino mass 
spectrum at the future T2K (Tokai-to-Kamiok a) |38l and 
NO^A J39j long-baseline experiments, see Ref. J40|. If the 
value of 9i3 turns out to be very small, the ultimate high- 
precision measurements may require super-beams |41|, 
neutrino factories 0> 03 i beta beams Q , or a combi- 
nation [i M Ii ^ . 

Despite the great promise of planned and future 
neutrino oscillation experiments, two essential neutrino 
properties would still remain unknown. Is the neutrino 
distinct from its antiparticlc (y ^ D, Dirac case), or iden- 
tical to it [y = D, Majorana case)? 9 What are the abso- 
lute values of the neutrino masses? 

Searches for neutrinoless double beta decay, a rare — 
and hitherto unobserved — transition between two nuclei 
with the same mass number (A) that changes the nuclear 
charge {Z) by two units, 

(Z 1 A)^(Z + 2,A) + e^ +e 2 , (15) 

are for now our only probe for Majorana neutrinos. Ob- 
servational upper limits on /3/3(Oi/) rates provide an upper 
bound on the so-called "effective Majorana mass" of the 
electron neutrino 

(m cff ) < 0.3- 1.0 eV . (16) 

Forthcoming f3[3{Qv) experiments that aim for sensitivity 
approaching (m e ff) ^ 0.05 eV could well establish neu- 
trinos as Majorana particles |5lll52l|. 

Direct information on the absolute scale of neutrino 
masses can be extracted from kinematical studies. The 



7 We consider here the canonical picture of three neutrino families. 
The ongoing MiniBooNE experiment l2ll is expected to explore 
all the parameter s pac e of the <-> i/ e mutation claimed by the 
LSND experiment kill , with its implication of additional (pre- 
sumably sterile) neutrinos. The neutrino mixing matrix is some- 
times called the Pontecorvo-Maki-Nakagawa— Sakata (PMNS) 
matrix. 

8 We define Am?. = m 2 — m 2 , where rrtj is the mass of vi. 



9 Prevailing theoretical opinion favors the Majorana character, if 
only because we know no princi ple t hat forbids Majorana mass 
terms. The see-saw mechanism l47l l48t 141.1 offers a natural in- 
terpretation of the smallness of neutrino masses, and points to a 
new scale associated with the heavy Majorana neutrino. The see- 
saw mechanism implies lepton-number-violating processes that 
could have happened in the early Universe. The decay of the 
heavy, weak-isoscalar, Majorana neutrinos, together with B + L- 
violating sphaleron transitions can give rise to the baryon asym- 
metry of the Universe loll . 
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FIG. 4: Composition of the neutrino mass eigenstates (^1,^2,^3) in terms of the flavor eigenstates {v e ,v,j,,v T ) for the normal 
(left panel) and inverted (right panel) hierarchies [from Ref. |3^]. Each bar shows the flavor mixture of a mass eigenstate as 
the CP- violating phase, 5, varies from cos 5 = — 1 (bottom) to cos 5 = 1 (top). The other mixing parameters are held fixed at 
the representative values sin 2 #12 = 0.30, sin 2 #13 = 0.03 and sin 2 823 = 0.50. To an excellent approximation, Am 2 2 = Am 2 ol 
and Am! 3 = Am 2 tm . 



present upper bound on the electron-neutrino mass from 
tritium beta-decay experiments is 2.2 eV (95% CL) [^3, 
l&jf . and in the future the KATRIN experiment is ex- 
pected to be sensitive to electron-neutrino masses ap- 
proaching 0.2 eV (90% CL) [55|. 

The evolution of the Universe is sensitive to the ab- 
solute neutrino mass scale, independent of mixing pa- 
rameters or CP-violating phases [561 l57| . for example, 
through the influence of neutrinos on the matter power 
spectrum jH^, Combining WMAP observations of 

the cosmic microwave background with large-scale struc- 
ture data from the 2dF Galaxy Redshift Survey [6(| or 
the Sloan Digital Sky Survey |6l| yields impressive con- 
straints on neutrino mass J2i m vi ~ 1 e V, the precise 
value depending on the priors and the data set 0, • 
However, the neutrino mass limits arising from existing 
large-scale structure can be evaded if new interactions 
(such as the coupling of neutrinos to a light boson) en- 
able the relic neutrinos to annihilate at late times 63] . 

In Figure [S] we show the allowed ranges of neutrino 
masses (mi,TO2, TO3) in terms of the lightest neutrino 
mass mi for the normal and inverted hierarchies, using 
the best- fit values of the solar and atmospheric mass gaps 
reported above. For the entire range of permitted masses, 
at least two relic species are nonrelativistic (m v ^> {p v o)) 
in the present Universe. Consult Ref. [65L l66| for re- 
cent surveys of the prospects for determining the absolute 
mass scale of the neutrinos. 



D. The Absorption Cross Section 



sources, the reactions of interest are usually the deeply 
inelastic scattering processes vN — > fi + anything or 
vN — > v + anything. Nucleons are so rare throughout the 
Universe at large (tib = (2.5 ± 0.1) x 10~ 7 cm~ 3 in the 
current Universe [3J) that neutrino-nucleon scattering is 
a negligible source of attenuation, even over cosmologi- 
cal distances. A path length of 8 x 10 5 Mpc in today's 
Universe corresponds to 1 cm water equivalent (cmwe). 
It is convenient to define the interaction length, 

= l/a vN {E v )N A , (17) 

where a v jss is the appropriate neutrino-nucleon cross sec- 
tion and N A = 6.02 2 x 10 23 moP 1 = 6.02 2 x 10 23 cm" 3 
(water equivalent) is Avogadro's number. For neutrino 
energies in the range 10 16 eV < E v < 10 21 eV, a recent 
calculation |63 yields (v, T>)N total cross sections 

/ E n0.363 

o-vjv—U w 0.78 x IO- 35 cm 2 (7-^) ■ (18) 

It is not unreasonable to extrapolate this form a few or- 
ders of magnitude higher in energy. 10 

We plot the resulting interaction length in Figure 
For the energies that will be of interest to us, the in- 
teraction length lies in the range 10 6 - 10 9 cmwe, or 
10 12 - 10 15 Mpc in the current Universe. These distances 
are extraordinarily vast, in view of the expectation that 
astrophysical sources of ultrahigh-energy neutrinos such 
as active galactic nuclei lie within 100 Mpc of Earth. At 



In terrestrial experiments that seek to detect neutri- 
nos originating in particle accelerators or astrophysical 



For an examination of different extrapolations in energy and of 
the influence of exotic mechanisms, see Ref. l68l . 
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FIG. 5: Favored values for the neutrino masses as functions of the lightest neutrino mass, me, in the three neutrino scenario 
for normal hierarchy (left panel, mi = mi) and the inverted hierarchy (right panel, me = 7713). [After Ref. |r34|.] 
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10 100 1000 10* 10° 10° 10' 10° 10" 10 lu 10" 10' 
E„ [GeV] 

FIG. 6: Interaction length defined in Eq. (11711 for the reac- 
tions vN — > anything as a function of the incident neutrino 
energy. The left-hand scale, in cmwe, is appropriate for ter- 
restrial applications; the right-hand scale, in Mpc for the cur- 
rent Universe, is appropriate for transport over cosmological 
distances. [After [6711 .1 



earlier — but not too early — epochs, the nucleon density 
scales with redshift as = n#(0)(l + z) 3 . Even 

so, back to z 20 when the first astrophysical neutrino 
sources began to shine, the interaction length is far too 
long for vN scattering to be a significant mechanism for 
reducing the flux of neutrinos incident on Earth. 11 



Over the energy range of interest for neutrino astron- 
omy, the interactions of v e ,v^,v T ,v^,v T with electrons 
can generally be neglected compared to interactions with 
nucleons. The case of v e e interactions is exceptional, be- 
cause of the intermediate-boson resonance formed in the 
neighborhood of E™ les = M^/2m e « 6.3 x 10 15 eV. 
The peak cross section, a(y e e — > W~ — > anything) ss 
5 x 1CP 31 cm 2 , corresponds to an interaction length at 
resonance of 6 x 10 6 cmwe [6^. Assuming that the den- 
sity of electrons throughout the current Universe is com- 
parable to the density of baryons, even resonant D e e scat- 
tering contributes a negligible attenuation of astrophysi- 
cal neutrinos en route to Earth. 

The number density of relic neutrinos (of each species) 
in the current Universe is some 2.2 million times the 
number density of baryons. The thicker relic neutrino 
target (!) combined with an appreciable cross section 
for vv — > Z° annihilation accounts for the importance of 
resonant absorption as an attenuator of extremely high 
energy neutrinos. As we will see shortly, the interaction 
length for resonant annihilation in the current Universe 
is C^ t 1.2 x 10 4 Mpc, some six orders of magnitude 
shorter than the interaction lengths for vN or resonant 
v e e scattering. 

The cross section for neutrino- antincutrino annihila- 
tion into fermion pairs through the Z° is given by 



s.Z 



(y a v a -> //) 



G 2 F miE v N t 



(/) 



3- 
L2 f 



(19) 



Rj 



[i-yy + Tl/Ml 



11 A complete treatment, taking into account redshifting of the neu- 



trino energy, does not change the conclusions |^ 
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TABLE I: Decay modes and branching fractions of the Z° Q. 



■30 



Decay mode 


Branching fraction 


hadrons (uu + dd + ss + cc + bb) 


70% 


charged leptons (e + e~ + fi + yT + t + t~ 


) 10% 


invisible (v a v e + v^v^ + v T u T ) 


20% 



where y = 2m i E 1/ /M%, G F = 1.166 37 x 10~ 5 GeV~ 2 is 

the Fermi constant, m, is the mass of the target (relic) 

( f) 

neutrino, E v is the incident neutrino energy, and Nc 
is the number of colors of the fcrmion /: 1 for leptons 
and 3 for quarks. The chiral couplings of / are Lf = 



.(f) 



2Qf sin 9w and Rf 



-2Q/siir 
If) . 



9\y, where Qf is 

the fermion's electric charge and t% ' = ±1 is the third 
component of its (left-handed) weak isospin. 12 

We have written Eq. (|19fl for the annihilation of a neu- 
trino of flavor a on its antincutrino counterpart. In our 
application, the cross section must be weighted by the 
probability for the mass eigenstate Vi to contain the fla- 
vor component v a , which is to say, by the absolute square 
of the appropriate neutrino mixing matrix element. We 
have also assumed that the neutrinos are Majorana parti- 
cles; the difference between Majorana and Dirac particles 
is explained in the following fll El 

When summed over the kinematically accessible de- 
cay products of Z , namely the charged leptons e, fi, r, 
the neutrinos v e ,v^,v T , and the quarks u,d,s,c,b } i|19|) 
leads to the branching fractions collected in Table [I] In 
our study of absorption lines, we shall regard the 20% of 
Z° — > vv decays as removed from the incident beam. In a 
detailed study of particular experimental circumstances, 
one might choose to improve this approximation. 

Depending on the incident beam and the relic tar- 
get, other processes may contribute. A complete cata- 
logue was given by Roulet |20|. whose notation we em- 
ulate here. We present the components of the neutrino- 
(anti)neutrino cross sections in Figure Neutrino- 
antincutrino scattering in general receives a contribution 
from t-channel Z-exchange, 



a ' (v a v p -> v a v p ) 



G\miE v 



where Fi(y) = [y 2 + 2y - 2(1 + y) ln(l + y)}/y 3 . The 



s-t interference term is 



G 2 F miE v 



2n 

(y-i) 



F 2 (y) (2i) 



{i-yy + vl/Ml 



12 We take M z = 91.1876 GeV, T z = 2.4952 GeV, sin 2 6 W = 
0.231; we have taken account of fermion masses in our numerical 
studies. 
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FIG. 7: Total neutrino annihilation cross section and the dif- 
ferent contibution channels as a function of the ultra- high neu- 
trino energy assuming a relic neutrino mass of ra v = 10 -5 eV 
and zero redshift. 



where F 2 (y) = [3^ 2 + 2y - 2(1 + yf in(i + y)]/y 3 . 

Neutrino-antineutrino scattering to a pair of charged 
leptons may proceed by W exchange in the t channel, 

<r tiW {u a Pp^£j fi ) = 4G2pmiEu F 1 (y). (22) 

For charged-lepton pair production, the interference be- 
tween the s-channel Z exchange and the i-channel W 
exchange is 

^G\m % E v 

. I, .1 I Wen.- - — , 



wz 



(y a v p -> l a l p ) = 5 a p- 



7T 

(y-i) 



-F 2 (^)(sin 2 Bw-h) 



(23) 



[l-yf + Ti/Mi 

Neutrino-neutrino (or antineutrino-neutrino) elastic 
scattering is mediated by i-channel Z exchange, with a 
cross section 

G\m i E v 1 



t:Z 



(v a V0 -* V a Vf3) = 



i + y 



(24) 



that is accompanied, for identical species, by the u- 
channcl contribution 



'{v a Vp- V a Vf3) = 5 a (3 



G 2 F m % E v 



y y{i + \y) 



(25) 



1 



Above the thresholds for W + W~ and Z°Z° pair pro- 
duction, we include the vv — > vector-boson-pair cross 
sections in our numerical analysis. The effect of these 
processes on neutrino attenuation is minor; the relevant 
formulas may be found in Ref. 20] . 
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E. Dirac versus Major ana Relics 

The interaction cross section may depend on whether 
the relic (target) neutrinos are Dirac or Majorana par- 
ticles. If the relic neutrinos are extremely relativistic, 
the Dirac and Majorana characters are indistinguish- 
able: relativistic neutrinos are pure left-handed chirality 
states, because only such states are produced in the weak 
interactions. Chirality and helicity coincide, and the 
right-handed chirality is absent. Nonrelativistic Dirac 
and Majorana neutrinos exhibit distinctive behavior. In 
the static limit, Dirac neutrinos are left-handed helic- 
ity eigenstates with equal populations of left- and right- 
handed chiralities; Dirac antineutrinos are right-handed 
helicity eigenstates, also with equal populations of lcft- 
and right-handed chiralities. Since only the left-handed 
neutrino (right-handed antineutrino) chiralities interact, 
the other two components of the Dirac neutrino field are 
sterile. Because Majorana neutrinos are their own an- 
tiparticlcs, both chiralities interact. Accordingly, in the 
static limit, the interaction cross section on a Majorana 
target is twice the cross section on a Dirac target. 

The interactions with a Majorana neutrino's "wrong- 
chirality" population enter with weight m 2 /(£„ + p„) 2 , 
where p v is the relic neutrino momentum and e„ = 
yjp' 2 , + to 2 ,, so that 



^Majorana 



= 1 



(26) 



CTDirac {Eu + Pv) 2 l+p v /e v 

The interaction length for annihilation on relic neutri- 
nos illustrates the difference between the Majorana and 
Dirac cases and the transition from static to extreme rel- 
ativistic regimes. We define 



• t -int:M,D 



(27) 



where E^ ICS = M^/2m u , and evaluate the ratio i|26|) 
using the mean momentum (p u ) and energy (e„) = 
({pD + rn 2 ,) 1 / 2 from the Fermi-Dirac distribution 0. 
In Figure [8] we depict the Majorana and Dirac interac- 
tion lengths over the range of redshifts considered in this 
study, for three illustrative values of the relic neutrino 
mass. 

Contrary to common wisdom, the distinction between 
Majorana and Dirac relics is not readily observable in 
neutrino absorption lines, unless the relic neutrino mass 
approaches « 0.1 eV, a value close to the present cos- 
mological upper bounds for a quasi-degenerate neutrino 
spectrum. Appreciably lighter relics are relativistic over 
much or all of the redshift range we consider. 



F. Super-High-Energy Neutrino Sources 

1. General orientation 



sub-eV relics that current information on the neutrino 
spectrum leads us to expect, the cosmic neutrinos must 
have energies no less than those of the highest-energy 
cosmic rays ever observed [z3,0|- Indeed, the required 
cosmic neutrinos have been named Super- GZK neutri- 
nos , as their energies lie above the Greisen-Zatsepin- 
Kuzmin |73l cutoff in the cosmic ray spectrum. 

It is worth taking a moment to review the GZK argu- 
ment, because it implies the existence of so-called cos- 
mogenic neutrinos |75j . Extremely high energy cosmic 
rays — let us take protons, to be concrete — can lose energy 
by interacting with the cosmic microwave background 
whose properties were recalled in fllBI The key energy- 
loss mechanism is pion photoproduction, p + 7 — > 7r + N. 
Ultrahigh-energy 1^,7, and arise from the decays of 

7T + , 7T°, 7T~. 

Taking the energy of a typical CMB photon as 



(Pio) 



-C(4) rp 

3 C(3)' T ° 



2.701T « 6.341 x 1CT 4 eV , (28) 



we estimate the threshold for pion production to be 



E i- 



run (m w + 2M p 
4(p 7 o) 



1.1 x 10 20 eV , 



(29) 



where is the pion mass and M p is the proton mass. 
Accordingly, any proton with energy > 10 20 eV that tra- 
verses a long path in the current Universe will suffer 
energy loss through pion photoproduction. The inter- 
action length, determined by scattering at the A(1232) 
resonance, is approximately 10 Mpc, short compared 
with the 100-Mpc distance to active galaxies. At ear- 
lier epochs, the interaction length (at redshifted energy) 
scales with the number density of CMB photons. 

No experiment has yet detected neutrinos with energies 
above 1 TeV that originate outside Earth's atmosphere. 
To discuss possible sources 13 of cosmic neutrinos that 
might be useful for absorption spectroscopy we enter a 
largely unexplored realm of upper limits and models not 
disciplined by extensive data sets. 

Both acceleration mechanisms and top-down (decay) 
phenomena may be at the origin of super-high-energy 
neutrinos. We consider these two classes of sources briefly 
in turn. 

Extragalactic objects such as active galactic nuclei 
(AGNs) and gamma-ray bursters (GRBs) are gener- 
ally regarded as promising sites for the production of 
ultrahigh-energy neutrinos. Protons accelerated to ex- 
treme energies may collide with the surrounding matter 
or the bath of photons to produce pions through the in- 
clusive reactions 



p + (N, 7) — » 7r + anything. 



(30) 



Observing absorption lines on the relic background re- 
quires an adequate neutrino flux at the resonant ener- 
gies E^ ICS = M%/2m Vi » 4.2 x 10 21 eV/m Vi . For the 



13 See the extensive review of cosmic-ray sources in Ref . l76l . 
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FIG. 8: Interaction lengths defined in Eq. 127H versus redshift at the Z resonance for neutrino masses m v = 10- 5 ,10- 3 ,10 _1 eV 
(left, center, and right panels). The left-hand scales are in centimeters, the right-hand scales in megaparsecs (1 Mpc = 
3.085678 x 10 24 cm). In the center and right panels, the lower (black) line is for the Dirac-neutrino case; the upper (red) line 
applies to Majorana neutrinos. 



If 7r + ,7r°,7r are produced in equal numbers, then the 
decay chains ir° — ► 77 and 

(and similarly for 7r~) imply products in the proportions 
7 : e + : e~ : : : v e : v e :: 2 : 1 : 1 : 2 : 2 : 1 : 1. If 
such processes were responsible for the flux of ultrahigh- 
energy gamma rays, then a similar flux of ultrahigh- 
energy neutrinos would be essentially guaranteed. 1 

The top-down scenarios — superheavy relic particles or 
topological defects formed in symmetry-breaking phase 
transitions predicted by unified theories, for example — 
do not require regions in which astrophysical processes 
can accelerate particles to super-high energies, but they 
depend on physics b eyon d the standard model that has 
not been established |8l|. They might populate energies 
beyond the reach of even the most extreme astrophysical 
processes, conceivably exceeding the scale on which the 
SU(3) c ®SU(2) L £g>U(l)y standard-model interactions are 
unified. 

In this class of models, ultrahigh-energy cosmic rays 
can be decay products of some supermassive X-particles 
with masses Mx close to the GUT scale. The supermas- 
sive X-particles could be long-lived relics of the early 
universe or could themselves arise from the collapse of 



A rather restrictive upper bound (Waxman-Bahcall) follows if 
energetic protons escape freely from suc h a source l77l . A more 
permissive (by ~ 40 x) cascade limit 173 relates the n eutr ino flux 
to the 7-ray flux observed by the EGRET instrument l79l aboard 
the Compton Gamma-Ray Observatory, assuming the photons 
are not obscured. Future 7-ray detectors, such as GLAST |8(_I . 
will improve the photon-flux baseline. In the case of hidden 
sources, from which neither nucleons nor photons escape, there 
is no way to bound the neutrino flux from above. 



topological defects. The X particles can decay into nucle- 
ons, gamma rays, and neutrinos with energies approach- 
ing Mx- 15 In the simplest models, the fluxes of neutri- 
nos that arise in this manner are bounded by the cas- 
cade limit. Hidden-sector topological defects, which may 
arise in multi-brane scenarios, evade the cascade limit 
and so might provide the largest (which is to say, least 
constrained) flux of Super-GZK neutrinos |8^| . 



2. Parametrizations of neutrino spectra 

Calculation of the cosmogenic neutrino flux is accom- 
plished by propagating an assumed primary proton flux 
through the cosmic microwave background over the rele- 
vant history of the universe, by means of transport codes. 
A standard Ansatz [82l l&jf is a power-law shape for the 
injection spectrum per unit comoving volume, 

<p p (E, z) = M ■ E- a (1 + z) m Q(E max - E) , (31) 

during the era characterized by z m { n < z < z max . Here 
B max represents the maximum energy to which protons 
can be accelerated by astrophysical processes, a is the 
spectral index, m is the redshift-evolution index, and M 
specifies the normalization. 

As an example, we show in Figure|Hl(£ , ^x) the cosmo- 
genic neutrino flux j(E u ) tuned [83 to saturate the cas- 
cade limit, j(E v )E% < 450 eV cm~ 2 s -1 sr -1 , derived 
from a recent analysis of the EGRET data [8(j. As an- 
ticipated, the cosmogenic neutrinos lie squarely in the 
domain of the GZK cutoff on the cosmic-ray spectrum. 
The flux depicted in Figure [5] is consistent with direct 
limits on the neutrino flux, which are summarized in 



See Ref. |82| and the works cited there for a general discussion. 
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FIG. 9: Solid line: Maximal cosmogenic neutrino fluxes per 
flavor computed for an injection spectrum characterized by 
£ max = 2 x 10 13 GeV, a = 1, m = 3, for < z < 2 (from 
Ref. H3). Dashed (red) line: Diffuse neutrino flux from hid- 
den topological defects in the form of necklaces, as calculated 
in Ref. |85| for a superheavy particle of mass M x = 10 14 GeV. 



Refs. |23j,|82(. Equal fluxes (at Earth) of all neutrino fla- 
vors is implied by the pattern of neutrino mixing, as we 
elaborate in the opening paragraphs of iflTl We comment 
briefly on event rates for planned detectors in illll CI 

To illustrate top-down scenarios, we cite a calculation 
of the diffuse neutrino flux arising from one species of 
topological defect, a necklace of monopoles and anti- 
monopoles strung on a cosmic string. Once formed in 
symmetry-breaking phase transitions in the early uni- 
verse, topological defects can survive indefinitely, un- 
less they collapse or annihilate |87j . producing massive 
quanta generically labelled X particles. As monopoles 
and antimonopolcs on a necklace annihilate, they pro- 
duce X particles in the form of heavy Higgs bosons and 
gauge bosons. These in turn are the source of ultrahigh- 
energy cosmic rays and neutrinos; analytic expressions 
for the fluxes in several cases are computed in Ref. [8l| . 

The dashed (red) curve in Figure shows the diffuse 
neutrino spectrum in the present universe that arises for 
monopole-antimonopole annihilations into particles with 
mass Mx — 10 14 GeV [85|. In top-down scenarios, ener- 
getic neutrinos might have been generated at very early 
times; in that case our computation of the absorption 
lines will entail an integration over redshift back to those 
early times. 



the detection of neutrinos from Supernova 1987A and 
from the Sun Jg^, . The detailed observation of at- 
mospheric neutrinos with GeV energies was crucial to 
establishing neutrino mixing and the AMANDA ex- 
periment has detected atmospheric neutrinos up to about 
10 5 GeV [gjj. Current exploration is dedicated to the 
search for extraterrestrial neutrinos — either from the dif- 
fuse background of active galactic nuclei or from point 
sources such as gamma-ray bursters — with energies be- 
tween 10 5 and 10 9 GeV. In addition to the ice-Cherenkov 
detector AMANDA II, dedicated neutrino telescopes in- 
clude the BAIKAL water- Cherenkov array [9l| and the 
antenna array RICE j^], which aims to detect ra- 
dio pulses emitted by neutrino-induced showers in the 
Antarctic ice. The Fly's Eye 93. HiRes [H, and 
AGASA [95| air-shower arrays are sensitive to horizontal 
air showers initiated by neutrino interactions deep in the 
atmosphere. 

The Goldstone Lunar Ultrahigh-energy neutrino Ex- 
periment (GLUE) [9(j has begun to search for radio 
emission from ultrahigh-energy cascades induced by neu- 
trinos or cosmic rays skimming the moon surface. The 
FORTE (Fast On-orbit Recording of Transient Events) 
satellite [97j has set upper limits on the UHE neutrino 
fluxes at energies beyond the GZK cutoff, looking for 
radio pulses generated by neutrino interactions in the 
Greenland ice sheet. 

Over the next few years, significantly increased sen- 
sitivities will be attained in IceCube [9^, a cubic- 
kilometer-scale ice-Cherenkov detector evolved from the 
AMANDA experience that is be ginn ing const ructi on at 
the South Pole. The ANTARES '99], NEMO [lO^, and 
NESTOR 101] projects are developing techniques for a 
cubic-kilome ter w ater-Cherenkov array in the Mediter- 
ranean Sea |102| . It is highly desirable that the next 
generation of neutrino telescopes not only characterize 
the incident neutrino energy and direction, but also tag 
the neutrino flavor by identifying the outgoing charged 
lepton. An optimistic assessment of prospects for flavor 
taggi ng i n detectors such as IceCube can be found in 
Ref. [io^ . 

ANITA |l04j . a balloon-borne array of radio anten- 
nas, will circle the Antarctic continent at an altitude of 
<~ 35 km to record radio bursts from neutrino interac- 
tions in the polar ice cap. The Pierre Auger Observa- 
tory, a 3000-km 2 hybrid detector for air showers in Ar- 
gentina's Mendoza province, will be sensitive to the inter- 
actions of > 10 9 -GeV neutrinos in the atmosphere |105| . 
Space-based instruments su ch a s the Extreme Universe 
Space Observatory (EUSO) Hoi a nd the Orbiting Wide- 
angle Light-collectors (OWL) |107| would have an energy 
threshold near 10 10 GeV. 



G. Detectors for Cosmic Neutrinos 

Neutrino astronomy is moving into a new, and much 
anticipated, era. The value of neutrino observatories has 
been clearly demonstrated in the MeV range through 



II. A TOY EXPERIMENT 

As a prelude to our investigation of neutrino absorp- 
tion spectroscopy in the physical Universe, we describe 
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a highly idealized situation in which an extremely high- 
energy neutrino beam traverses a very long column with 
the relic-neutrino properties of the current Universe. We 
neglect for now the expansion of the Universe and the 
thermal motion of the relic neutrinos. The "cosmic neu- 
trino attenuator" is thus a column of length L with uni- 
form neutrino density n„o = 56 cm -3 of each neutrino 
species, v e , v e , , , v T , v T . 

We imagine a terrestrial detector capable of dis- 
tinguishing arriving neutrino species through charged- 
current interactions producing electrons, muons, or tau 
lcptons, and of inferring the energies of the neutrinos that 
initiated those interactions. A detector of the required 
scale is unlikely to measure the charge of the outgoing 
lepton, and so would not distinguish neutrinos from an- 
tineutrinos in the beam. Our aim here is to identify the 
sorts of observations that might be made, should neu- 
trino absorption line spectroscopy become practical, and 
to point out how various neutrino properties would man- 
ifest themselves. 

We assume that the incoming neutrino beam orig- 
inates at least 100 Mpc from Earth, that it contains 
z/ e , P e , , P„, z/ T , v T in sufficient numbers to allow a mea- 
surement of the energy spectra of the neutrinos arriv- 
ing at Earth, and that the neutrino energy spectrum 
at the source is reasonably smooth. 16 Neutrino oscil- 
lations tend to produce a "beam" with roughly equal 
mixtures of u e ,u^,u T , whatever the flavor mix ture of th e 
extremely high-energy neutrinos at the source |l08l Il09| . 
The vacuum oscillation length, L osc = 4irE,s/\Am 2 \, 
is typically short compared with the intergalactic dis- 
tances we contemplate. For |Am 2 | = Am 2 2 as 8.2 x 
10 -5 eV 2 , for example, the oscillation length is L osc w 
3 x 10~ 25 Mpc • [E v /1 eV). At the resonance energy 
£jf res = Mz/2m v , the (solar) oscillation length is thus 
L osc w 1.2 x 10~ 3 Mpc/(m v< /l eV). Only for the nor- 
mal hierarchy with mi < 10~ 4 eV is the oscillation 
length not a negligible distance. The decoherence length 
for neutrinos above 10 18 eV is many orders of magni- 
tude greater than the Hubble distance, Dh = c/Hq = 
4200 Mpc(0.7//i) 0. 



A. An Idealized Experiment 

As we reviewed in ill CI the mass mi of the lightest neu- 
trino, which corresponds to m\ for the normal spectrum 
or to 7713 for the inverted spectrum, lies between 10 -5 eV 
and 0.1 eV. We present here studies for three example 
values: mi = 10~ 5 , 10~ 3 , and 0.1 eV. For a given mass 
of the lightest neutrino, the graphs in Figure \5\ show our 
current expectations for the masses of the other neutri- 
nos. 



If the column of relic neutrinos is thick enough to at- 
tenuate neutrinos appreciably through resonant absorp- 
tion at the Z° gauge boson, the energies that display 
absorption dips point to the neutrino masses through 
the condition m v = M^/2E^ res . 17 Evidently the heavi- 
est neutrino corresponds to the lowest-energy absorption 
line. For a normal hierarchy, we expect three or two or 
one absorption lines, as mi increases; for an inverted hi- 
erarchy, two or one absorption lines. 

We illustrate this prospect in Figure ^| where we plot 
the survival probabilities for v e ^v^, v T over the relevant 
range of neutrino energies, for me = 10~ 5 and 10 -3 eV. 
We imagine that the survival probability can be esti- 
mated reliably by fitting the shape of the energy spec- 
trum away from the dips, in much the same spirit as 
backgrounds are estimated in experiments that observe 
peaks. The amount of attenuation is governed by the 
length of the column and by the flavor composition of 
the neutrino mass eigenstates, which was indicated in 
Figure 0] 

As in other aspects of ultrahigh-cnergy neutrino stud- 
ies, flavor tagging at the d etector greatly enriches the 
scientific program ESlliiil. With flavor tag gmg, sim- 
ply looking for variations in the relative fluxes v e /al\, 
Vfj,/aSl, /all with energy may be highly revealing. We 
plot in Figure ITT1 the flux ratios v e jv^ and v T /v^ for nor- 
mal and inverted hierarchies with mi — 10~ 5 eV. 18 The 
flavor ratios are a powerful discriminant between the nor- 
mal and inverted hierarchies, because the v^jv^ ratio is 
an excellent diagnostic for the mass eigenstates — even if 
only the lowest-energy (highest-mass) dip is visible. In 
the normal hierarchy, muon neutrinos are depleted with 
respect to electron neutrinos, while for the inverted hier- 
archy, v e is depleted with respect to v^. This distinction 
can be essential in case the incident neutrino flux runs 
out before all of the absorption dips are revealed. 

If the neutrino masses exhibit a distinct hierarchy — 
as is the case for mi < 10~ 3 eV — then it is possi- 
ble, with adequate flux, to distinguish the normal and 
inverted hierarchies, even without flavor identification. 
In Figure 1121 we show the all-flavor survival probabil- 
ities, for a representative relic neutrino column length 
of 10 5 Mpc, for the normal and inverted spectra with 
mi = 10~ 5 , 10~ 3 , 10" 1 eV. For the normal hierarchy, we 
observe three, three, and one absorption lines; for the 
inverted hierarchy, two, two, and one dips. 



B. Majorana or Dirac Relics? 

We reviewed in ill El the implications of Majorana ver- 
sus Dirac character for the annihilation cross sections, 



It is a reasonable bet, though not essential to our analysis, that 
the beam contains an equal mix of neutrinos and antincutrinos. 



See Ref. IllOl for a related strategy to determine the neutrino 
masses. 

We assume the natural mix, v e : : u T :: 1 : 1 : 1. 
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FIG. 10: Survival probabilities for v e , v^, and v T in the ideal experimental scenario, for a normal hierarchy with me — 10 5 eV 
(left panel) or me = 10 -3 eV (right panel). The red and blue curves correspond to L = 10 4 and 10 5 Mpc. 
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FIG. 11: Flux ratios Ve/Vy, and Ur/Vfj, at Earth, for normal (left panel) and inverted (right panel) mass hierarchies with 
me = 10" 5 eV. The column length is L = 10 5 Mpc. 



and showed in Figure [S] how the distinction depends on 
neutrino mass and temperature. Absent a calibration of 
the length of the relic-neutrino column, we do not see how 
to exploit the factor-of-two difference between Dirac and 
Majorana cross sections in the static limit to distinguish 
the two cases. Even in the idealized (zero neutrino tem- 



perature) situation we treat in this Section, it is hard to 
imagine inferring the neutrino character from the depth 
of a single absorption line. If the lightest neutrino were 
relativistic and the heaviest neutrino nonrelativistic, it 
is remotely possible to imagine calibrating the effective 
column density on the highest-energy (lowest neutrino 
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mass) dip, and then noticing a smaller apparent depth 
for the lowest-energy (highest neutrino mass) dip, in the 
case of Dirac relics. In the more realistic circumstances 
we shall describe below, with integration over rcdshift 
and attention to the thermal motion of the relics, we do 
not see how to distinguish Dirac from Majorana relics. 



C. Unstable Relics 

To this point, we have considered neutrinos as sta- 
ble particles. "Invisible" decays, such as the decay of 
a heavy neutrino into a lighter neutrino plus a very 
light — or massl ess-(pseudo)scalar boson such as the ma- 
joron are not very well constrained by observa- 

tions |103]. 19 If CPT invariance holds, SN1987a data set 
an upper limit on the lifetime of the lightest neutrino of 
Tg/mg > 10 5 s/eV. Observations of solar neutrinos lead 
to T2/TO2 ^ 10~ 4 s/eV. Finally, if the neutrino spectrum 
is normal, the data on atmospheric neutrinos coming up- 
ward through the Earth, imply T^jra-j, > 10~ 10 s/cV. All 
of these bounds leave open the possibility that heavy relic 
neutrinos might have long since decayed away, so that the 



current Universe is filled not with the primordial neutri- 
nos, but with their decay products. 20 

Let us consider, for illustration, the simplest case of un- 
stable neutrinos, in which the two heavier neutrinos have 
decayed into the lightest neutrino plus invisible products. 
In that event, the density of the heavier neutrinos is now 
zero, but the density of the lightest neutrino is three 
times n„o- We further assume that the neutrinos in our 
beam have not decayed in flight. We show in Figure IT51 
the neutrino survival probability for the idealized experi- 
ment, for the normal hierarchy with mg = mi — 10 -3 eV. 
In contrast to what we see in the right-hand panel of Fig- 
ure llUI there is only one absorption dip located at the res- 
onant energy that corresponds to the lightest neutrino — 
the surviving relic. The dip is more pronounced for fixed 
column density than in the standard picture, reflecting 
the threefold increase in the population of vg. 21 

If only the lightest neutrino survives as a relic, the ab- 
sorption line is necessarily at a higher energy than in the 
standard case of three surviving flavors. In the extreme 
case that mg = 10~ 5 eV, the single absorption line would 
occur at E v w 10 18 GeV, so the demands for adequate 
flux at the highest energies are very great. 



A majoron too massive to serve as a neutrino decay product can 
nevertheless have important consequences for cosmology, includ- 
ing deviations from the standard expectations for the radiation 
energy density and changes in the positions of pe aks in the cos- 
mic microwave background power spectrum 11141 . 



20 As we noted in ill CI majoron interactions that would mediate 
neutrino decays can also mediate (co)annihilations that might 
lead to a vanishing relic neutrino density today. 

21 We suppose that all the vi and V3 relics have decayed to v\, and 
neglect cosmic expansion. 
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FIG. 13: Survival probabilities for v a ,Vy,, and v r in the ideal- 
ized experimental scenario, for the normal hierarchy in which 
all the heavy relics (i>2 and ^3) have decayed. The mass of 
the surviving (lightest) relic is mi = mi = 1CP 3 eV. The 
curves correspond to column lengths L = 10 4 Mpc (red) and 
10 5 Mpc (blue). 



We note that, with flavor identification, even the single 
absorption line could signal the normal or inverted hier- 
archy. In the normal hierarchy, the v e flux is strongly 
absorbed, whereas in the inverted hierarchy it is atten- 
uated only slightly. Without flavor tagging, there is no 
prospect of unmasking the hierarchy, if the heavy neutri- 
nos are absent as relics. 



III. NEUTRINO ABSORPTION LINES IN AN 
EVOLVING UNIVERSE 

Even if we had the means to prepare neutrino beams of 
the requisite energy, the Universe would not hold still for 
us to perform the idealized measurements described in 
flTTl The time required to traverse one interaction length 
for vv — > Z° annihilation on the relic background in the 
current Universe (1.2 x 10 4 Mpc = 39 Gly) exceeds the 
age of the Universe, 22 not to mention the human at- 
tention span. If we are ever to detect the attenuation 
of neutrinos on the relic-neutrino background, we shall 
have to make use of astrophysical or cosmological neu- 
trinos sources traversing the Universe over cosmic time 
scales. The expansion of the Universe over the propaga- 



tion time of the neutrinos entails two important effects: 
the evolution of relic- neutrino density and the redshift of 
the incident neutrino energy. We consider both effects in 
this Section. 



A. The Influence of Expansion and Redshift 

We displayed the redshift dependence of the relic- 
neutrino number density in Figure |21 in fllBI There we 
also introduced the column density defined in Eq. (|14|) 
and plotted in Figure|3 which characterizes the neutrino 
number density per unit of redshift. 

As we look back in time, the present energy E v q of 
an ultrahigh-energy neutrino increases with redshift as 
E v {z) = (1 + z)E v q. This scaling is easily understood: 
the energy of an ultrarelativistic particle is inversely pro- 
portional to its wavelength, which stretches out in an 
expanding Universe as (1 + z) _1 . 

We evaluate the survival probabilities as 



V(E u0 ) = exp 



dza uv {{\ + z)E u0 ) 



dn Vi (z) 



dz 



, (32) 



where a vv is the interaction cross section of gTDjweighted 
with the appropriate mass-flavor mixing factor and 
dn u Jdz is the column density, Eq. (|14fl . 23 The atten- 
uation process is thus summed over redshift, with the 
simplifying assumption that the neutrino sources in ques- 
tion switched on at a time in the past characterized by a 
single redshift, which we take to be z v < 20. 

Plausible acceleration sources of ultrahigh-energy neu- 
trinos, including active galactic nuclei and gamma-ray 
bursters, lie at redshifts of a few [s3| - Extremely ener- 
getic neutrinos might have been generated in the decays 
of su per-heavy particles or topological defects |8lL l82l 183. 
1115j long before the first stars began to shine. Their in- 
teraction with the relic neutrino background could reach 
as far back as the epoch of light neutrino decoupling, 

z~o(io 10 ) m. 

We show in Figure ^] the survival probabilities for 
v e ,v^,i' T integrated from the present back to redshift 
z = 20, for mi — 10 -5 and 10 -3 eV. Comparing with 
the analogous idealized "toy experiment" whose outcome 
was depicted in Figure we find that the absorption 
dips are distorted by the redshift of neutrino energies. 
Neutrino-antincutrino annihilation into Z° at redshift z 
implies a depletion of neutrino flux now at an energy 



E, 



Mf 



2m Vz {l + z) 



(33) 



The correspondence between neutrino mass and energy 
of an absorption line is compromised. No longer does the 



Some 13 Gly according to the current best estimates. 



23 We choose not to smear the cross section over a neutrino spec- 
trum, to make our analysis as generally informative as possible. 
Equation 1321 would result for any power-law neutrino spectrum. 
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FIG. 14: Survival probabilities for z/ e , and v T after an integration back to redshift z = 20, for a normal hierarchy with 
mi = 10 -5 eV (left panel) or mi = 10 -3 eV (right panel), as a function of the energy of the UHE neutrino. The scale at the 
top shows the neutrino mass defined as m„ = M%/2E V that would be inferred if neutrino energies were not redshifted. 



un-redshifted relation m u = M^/2E^ res allow the neu- 
trino mass to be inferred reliably. In these examples, that 
simple relation would overestimate neutrino masses by 
about an order of magnitude. In compensation, moving 
the dips to lower energies may be an important advantage 
for an experiment that is almost sure to be flux-limited. 

If terrestrial experiments determine the neutrino 
masses — or if other considerations place increasingly 
stringent bounds on the neutrino masses — then the loca- 
tion of the neutrino absorption lines might provide new 
information about the neutrino sources. If, for example, 
the location of the dip lies much more than one order 
of magnitude below the anticipated energy, that would 
be strong evidence that sources of extremely high energy 
neutrinos existed before the stars — a surprising and im- 
portant conclusion. 

The morphology of the distortion is illustrated in Fig- 
ure ^| where we plot the v e survival probability binned 
in redshift for the simple (one-dip) case of a degenerate 
neutrino spectrum with mi — 0.1 eV. Evidently binning 
in redshift is not an observational possibility, but we can 
impose it on our calculation to deconstruct the origin 
of the shift and broadening of the absorption line. As 
expected, we observe that the downward shift in the at- 
tenuated energy in the present Universe grows with the 
redshift at which the annihilation occurred. Moreover, 
the depth of the absorption dips is greater at higher red- 
shifts, because of the increased column density of relic 
neutrinos, Eq. JT3J. 

In an evolving Universe, the neutrino flux ratios remain 
an effective discriminant of the mass hierarchy. Figure H^l 
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FIG. 15: Electron-neutrino survival probabilities integrated 
back to 2 = 20 for a degenerate neutrino spectrum with mi = 
0.1 eV. The thick black curve shows the distorted absorption 
line after integration back to redshift z — 20, and the colored 
lines (from right to left) show the contributions from bins in 
redshift: red, < z < 4; green, 4 < z < 8; blue, 8 < z < 12; 
yellow, 12 < z < 16; magenta, 16 < z < 20. 



shows the survival probabilities for an inverted spectrum 
with mi — 10~ 5 eV. We plot in Figure El the flux ra- 
tios f e /f M and Vt/vh for normal and inverted hierarchies 
with mi = 10~ 5 eV, integrated back from the present 
to redshift z = 20. Compared with the toy- model ra- 
tios plotted in Figure ^2 the features are broadened and 
displaced toward lower energies, but the essential mes- 
sage is unchanged. In the normal hierarchy, f e /V M > 1 at 
the lowest-energy line, whereas in the inverted hierarchy 
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FIG. 16: Survival probabilities for v e , and v T after an inte- 
gration back to redshift z = 20, for an inverted hierarchy with 
me = 10 -5 eV, as a function of the energy of the UHE neu- 
trino. The scale at the top shows the neutrino mass defined 
as m„ = Mz/2E V that would be inferred if neutrino energies 
were not redshifted. 



V e /Vu, < 1. 



Alternative Thermal Histories 



Since the redshift dependence of the relic-neutrino col- 
umn density is imprinted on the neutrino absorption 
lines, it might be possible to infer some information about 
the column density from the absorption lineshape. As we 
saw in the discussion surrounding Eqs. (|T2*|l - (fTljl . each 
particular cosmology influences the relic-neutrino column 
density through the Hubble parameter (Tffi . 

We have adopted the ACDM model as our standard 
reference cosmology. Let us compare its implications for 
neutrino absorption lines with those of its predecessor, 
the SCDM (for standard cold dark matter) model charac- 
terized by matter density f2 m = 1 and zero cosmological 
constant. 

In Figure 1181 we compare the relic neutrino column 
densities for these two cosmological models. The column 
density is sytematically larger in the ACDM model, be- 
cause flat-Universe models with a cosmological constant 
imply a larger physical volume associated with unit red- 
shift than flat-Universe models dominated by matter. As 
Figure H5l shows, neutrino attenuation is more efficient in 
the ACDM model than in the SCDM model. There, as in 
Figure 1151 we show our expectations for the v e survival 
probabilities integrated back to z = 20, for a degenerate 



neutrino spectrum with mi = 0.1 eV. 

If we had a "standard candle" for neutrinos analogous 
to the Type la Supernovae for photons, we might imagine 
discriminating between these two cosmological models, 
because the difference we compute here is similar to what 
is revealed by the supernova luminosity distances. Just 
as distant supernovae appear fainter in a dark-energy 
dominated Universe, so are the extremely high energy 
neutrinos more thoroughly absorbed. 

The discrimination among cosmological models would 
be more acute if extremely high energy neutrino 
creation — and attenuation — were initiated at much ear- 
lier times, or much larger redshifts, than we consider 
in this example. If supermassive particle decays are an 
additional — early — neutrino source, then different ther- 
mal histories might yield dramatically different neutrino 
absorption lines. 



C. Event rates 

Since experiments have just begun to explore the spec- 
trum of ultrahigh-energy neutrinos, quantifying the chal- 
lenge of establishing neutrino absorption lines is a very 
uncertain undertaking. We refer to §111 of Ref. j23| for a 
useful assessment. We concur that observing the cosmic- 
neutrino absorption lines with planned detectors will re- 
quire extended exposures, exotic sources of ultrahigh- 
energy neutrinos, or both. Here we briefly examine the 
case of quasidegenerate neutrino masses, which mini- 
mizes the energy required of the incoming neutrinos. 

Neutrino absorption lines in the flux of cosmo- 
genic neutrinos will appear at neutrino energies E v > 
10 12 GeV, whereas those in the flux of neutrinos gener- 
ated much earlier from topological defects may be red- 
shifted to energies E v < 10 10 GeV. Such energies lie 
in t he do main of radio- Cherenkov and air-shower detec- 
tors |l!6| . IceCube is optimized for neutrinos in the TeV- 
to-PeV range; it is sensitive to neutrinos of much higher 
energies, though the ability to characterize the neutrino 
energy di minish es progressively. IceCube and possible 
extensions 117] hold promise for identifying the flavor of 
the ultrahigh-energy neutrinos. We explore here the ex- 
pected sensitivit ies t o the absorption lines in the future 
ANITA mission [i"fJ4| . 

Suppose now that super-high-cnergy neutrinos arise 
from the decay of superheavy particles created in the 
collapse (or annihilation) of topological defects. If 
these neutrinos are produced at very early times, de- 
tecting their annihilation on the relic-neutrino back- 
ground may give us a rare glimpse of processes that 
prevailed in the very early universe. For exam- 
ple, the monopole-antimonopole boundstates dubbed 
"monopolonium" would have formed at the time of he- 
lium synthesis Neutrinos released in the decay of 
X particles, as described in ^1 F 21 have traversed many 
different epochs; the expansion of the universe during the 
neutrino's propagation will be imprinted on the neutrino 
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FIG. 17: Flux ratios Ve/v^ and Vrjv^ at Earth, for normal (left panel) and inverted (right panel) mass hierarchies with 
me — 10 _J eV, after integration back to redshift z = 20. The scale at the top shows the neutrino mass defined as m v 
that would be inferred if neutrino energies were not redshifted. 
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FIG. 18: Relic neutrino column densities versus redshift for 
the ACDM model (upper black curve) and SCDM model 
(lower red curve). 



spectrum. It remains to be seen whether we can learn to 
read the imprints in the absorption lines. 

We show in Figure [23 the survival probability for 
neutrinos emitted in the decay of 10 14 -GeV X particles, 
following Ref. |8l|. For a degenerate neutrino spectrum, 
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FIG. 19: Electron-neutrino survival probabilities integrated 
back to 2 = 20 for a degenerate neutrino spectrum with me = 
0.1 eV in the ACDM model (lower black curve) and SCDM 
model (upper red curve). 



the absorption line occurs at E v % 10 10 GeV, well within 
the range to be covered by future experiments, includ- 
ing ANITA and EUSO. In the most optimistic scenario 
for the TD neutrino fluxes, a 5-er significance level (164 
events) co uld b e achieved after a single fifteen-day flight 
of ANITA |104| . For these redshifted, early-time neutri- 
nos, we must face an additional complication: the pres- 
ence of relatively late-time neutrino sources, such as the 
cosmogenic mechanism, that provide a pedestal under 
the signal. If the foreground sources are too bright, com- 
pared to the flux from top-down sources, more events will 
be needed to establish the absorption lines. 
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FIG. 20: Survival probability for electron neutrinos created in 
the decay of a superheavy particle, X, with Mx = 10 14 GeV, 
for the case of mi = 0.1 eV. 



IV. EFFECT OF THE RELIC-NEUTRINO 
TEMPERATURE 

Relic neutrinos are moving targets, with their momen- 
tum distribution characterized in the present Universe by 
Eq. ©. The thermal motion of the neutrinos gives rise 
to a Fermi (momentum) smearing of the UHE-j/-relic-^ 
cross section. The resonant incident-neutrino energy for 
a relic neutrino in motion is given by 



2(£V -p v cost 



(34) 



where p v and e„ are the relic-neutrino momentum and 
energy. The angle 9 characterizes the direction of the relic 
neutrino with respect to the line of flight of the incident 
UHE neutrino. Accordingly, the resonant energy will be 
displaced downward from M^/2m v to approximately 



Mf 



(35) 



where {e v ) = [(p 2 ,) + m 2 ] 1 ' 2 plays the role of an effective 
relic-neutrino mass. We plot the effective mass in Fig- 
ure 1211 for the interesting range of neutrino masses and 
redshifts. 

The root-mean-square relic-neutrino momentum, 
which ranges from 6 x 10~ 4 eV in the present Universe 
to 2.5 x 10~ 2 eV at redshift z — 20, thus serves as 
a rough lower bound on the effective neutrino mass. 
At a given redshift, the resonance peak for scattering 
from any neutrino with m v < will be changed 

significantly. 

We display the effect of Fermi motion on the in- 
formation cross section in Figure 1221 The annihilation 
cross section depends on the relic neutrino energy and 
momentum. We have integrated the exact expression 
for the cross section weighted with the Fermi-Dirac mo- 
mentum distribution of the relic neutrinos over the relic 
neutrino momentum phase space. We show two series 
corresponding to redshifts z — and 20 and relic neu- 
trino masses 10~ 5 , 10~ 4 , 10~ 3 , 10~ 2 , and 10" 1 eV. In 
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FIG. 21: Effective masses of relic neutrinos, estimated by 
Ev — [ip^+m 2 ] 1 ^ 2 , as functions of the neutrino mass m„. The 
mean-squared neutrino momentum in the current Universe is 
given by Eq. 1101 . From bottom to top, the curves refer to 
redshifts z from to 20, in steps of 2. 



each panel, the narrow peak (red curve) shows the anni- 
hilation cross section as a function of incident neutrino 
energy for a stationary relic- neutrino target; the broad 
peak (blue curve) shows the annihilation cross section 
averaged over the thermal distribution of relic-neutrino 
momenta. In every case, the thermally averaged cross 
section peaks near E^f ros , indicated by a downward arrow. 
Consequently, the neutrino mass inferred from the loca- 
tion of an absorption line in a hypothetical experiment 
at fixed redshift would be approximately (e v {z)), rather 
than the true neutrino mass, whenever (e v (z)) ^ m v . 

Now we consider the influence of thermal motion on 
the structure of neutrino absorption lines in an evolving 
Universe, integrating over redshift as we did in ^III Al We 
show in Figure 1231 the survival probabilities for v e ,i/^, v r 
integrated from the present back to redshift z = 20, for 
mg = 10~ 5 and 10~ 3 eV. Comparing with the analogous 
calculation neglecting the Fermi motion of the relic neu- 
trinos, whose outcome was depicted in Figure El we nn d 
that the three distinct dips have been merged into a single 
complex dip. As we could anticipate from the displace- 
ment of the resonant energies displayed in Figure l2*2l the 
prospect of determining distinct absolute neutrino masses 
has faded appreciably. 

The distinction between the normal and inverted hier- 
archy remains, however. We plot in Figurc l2*4"l the survival 
probabilities in the case of an inverted hierarchy with 
me = 10 -5 eV. Again, Fermi motion has merged the two 
distinct absorption lines (compare Figure into one. 
The pattern of v e , Vu attenuation is different from the 
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FIG. 22: Neutrino-antineutrino (annihilation) cross section as a function of the incident neutrino energy. In each panel, the narrow peak (red curve) applies for a 
relic neutrino at rest. The broader (blue) curves, shifted to lower energies result when the Fermi motion of the thermal distribution of relic neutrinos is taken into 
account. The upper series corresponds to relics in the present Universe (z = 0); the lower series corresponds to redshift z — 20. From left to right, the panels depict 
relic neutrino masses 10~ 5 , 10~ 4 , 10~ 3 , 10~ 2 , and 10" 1 eV. Arrows designate the thermally displaced resonant energy E^ TeB given by Eq. 1351 . 
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FIG. 23: Survival probabilities for u e , v^, and v T as a function of the neutrino energy, after integration back to redshift z = 20, 
taking into account the Fermi smearing induced by the thermal motion of the relic neutrinos. The results apply for a normal 
hierarchy with lightest neutrino mass me = 10 -5 eV (left panel) or me = 10 -3 eV (right panel). 



normal-hierarchy case, as we show in Figure l2~5l In the 
normal- hierarchy case, v^Jv^ > 1 for the prominent dip 
at the lowest energy, whereas for the inverted hierarchy 

Ve/Vp < 1. 

Although the thermal motion of the relics is largely 
responsible for obliterating the three-dip structure, we 
might inquire whether some of the structure would sur- 
vive if the sources of ultrahigh-energy neutrinos came 
into existence more recently than redshift z = 20. In 
Figure ED we show the survival probabilities for normal 
and inverted hierarchies, with the lightest neutrino mass 
rri£ = 10 -5 eV, integrating from the present back to red- 
shift z = 10. The multiple dip structure is somewhat 
more pronounced than in the z — 20 scenario, but still 
indistinct. Because the Universe is more transparent 
to neutrinos when evolved back through z = 10 than 
evolved back through z = 20, the dips are less prominent 
than in Figures |2"S1 and |2"H Flux ratios again distinguish 
the normal and inverted hierarchies. 



V. UNCONVENTIONAL NEUTRINO 
HISTORIES 

A. Lepton asymmetry in the early Universe 

If the number of leptons differed from the number 
of antileptons in the early universe, the present den- 
sity of relic neutrinos might be very different from the 
expected value of 56 cm -3 . In the standard big-bang 
model, the neutrino asymmetry r\ Va = (n Ua — n Pa )/nj, 



is assumed to be negligibly small, comparable to the 
tiny-but-crucial baryo n asy metry t)b = (^s — n B)/ n -y = 
(6.1 ± 0.2) x 10~ 10 jUnf. 24 Such tracking is typical of 
scenarios in which B — L is a conserved quantum number 
and leptogenesis is the forerunner of baryogenesis. How- 
ever, several models are k nown in which a small -qs is 
accompanied by a large r\ v |ll9t ll20]. Consequently, it is 
worth our while to examine the observational constraints 
on the neutrino asymmetry. 

A significant lepton asymmetry in the early universe, 
represented by a nonvanishing chemical potential £ Q , 
would increase the expansion rate of the universe, in- 
creasing the amplitude of the acoustic peaks observed 
in the cosmic microwave background anisotropy. A 
combined analysis of CMB and SN la data |l2l| finds 
-0.01 < £ e < 0.22 and |^ )T | < 2.6 at 95% CL. (The v e 
asymmetry is more tightly constrained because it influ- 
ences the neutron-to-proton ratio at the epoch of big- 
bang nucleosynthesis, and hence the helium fraction.) 
If neutrino mixing equilibrates all the chemical poten- 
tials be fore BBN, then the bo und for £ R applies to all 
species |l22j . and |f„ a | < 0.1 |l23t Il24| . whereupon in 
the current universe n„ a k, 56 cm -3 . However, if a ma- 
joron field shoul d block n e utrin o oscillations in the pri- 
mordial plasma |l25l Il26t Il27| , then ^ and f T do not 
equilibrate with £ e . If we take the largest allowed val- 



The apparent electrical neutrality of the universe implies a sim- 
ilarly small electron-positron asymmetry. 
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FIG. 24: Survival probabilities for v e , v^, and v T as a function 
of the neutrino energy, after integration back to redshift z = 
20, taking into account the Fermi smearing induced by the 
thermal motion of the relic neutrinos. The results apply for an 
inverted hierarchy with lightest neutrino mass raj = 10 eV. 



ues, |£ M , r | ~ 3, then a large le pton asymmetry appears, 
n Va — n Pa ~ n Ua ~ 1050 cm -3 |l2Sj . The increased neu- 
trino density would reduce the interaction length for an- 
nihilation on the favored relic neutrino species. 



B. Neutrinoless Universe 

It is conceivable that neutrinos may interact in ways 
not foreseen in the standard electroweak theory. If, for 
example, the neutrinos couple to an extremely light bo- 
son, relic neutrinos might annihilate at late times, greatly 
reducing — indeed, likely erasing — the neutrino density 
expected in the standard cosmology with standard-model 
interactions |63j |. If no neutrino relics remain in the re- 
cent (z < 20) Universe, then no neutrino absorption lines 
will be observed. To demonstrate that the search for ab- 
sorption lines has covered the right terrain, it will be ex- 
tremely helpful to establish the absolute scale of neutrino 
masses. 



C. Neutrino Clustering 

One early inspiration for the study of ultrahigh-energy 
neutrino absorption was the prospect that massive neu- 
trinos, with masses in the range 30-100 eV, might ac- 
count for the dark matter in the Universe [2(i|. Had 
that been the case, neutrinos might have been absorbed 



not only on relics at large redshift, but also on relics 
clustered in galactic halos. For neutrino masses in the 
range we now consider likely, m„ < 0.1 eV, gravita tiona l 
clustering of neutrinos appears to be insignificant |129| . 
Moreover, since light neutrinos would cluster — or accrete 
onto other structures — only at late times, absorption 
spectroscopy is unlikely to record neutrino density con- 
trasts E3 ■ The question of detecting neutrino clumps by 
other techniques remains open. 



VI. SUMMATION AND OUTLOOK 

Let us assess the prospects for developing cosmic- 
neutrino absorption spectroscopy as a new window on 
Nature. First, we require that neutrino observatories es- 
tablish the existence of cosmic-neutrino fluxes that ex- 
tend to energies in the Greisen-Zatsepin-Kuzmin regime 
and beyond. If the requisite neutrino fluxes exist, then 
the essential task will be to observe absorption lines at 
energies corresponding to the masses of one or more neu- 
trino species. Making those observations will require 
detectors with vast effective volumes, or long exposure 
times, or both. Detecting the relic-neutrino background 
would, by itself, be an enormously satisfying observa- 
tional accomplishment. A survey of other tech niques that 
might detect the relics is given in Ref. [l28j. Short of 
direct detection, it may be possible to quantify the in- 
fluence of primordial anisotropics in the relic neu trinos 
upon the cosmic microwave background 

EMEU. 

The observation of cosmic-neutrino absorption lines 
will open the way — at least in principle — to new insights 
about neutrino properties and the thermal history of 
the universe. Our calculations, with their successive in- 
clusion of potentially significant effects, show that how 
the tale unfolds will depend on factors we cannot fore- 
see. The earlier in redshift the relevant cosmic-neutrino 
sources appear, the lower the present-day energy of the 
absorption lines and the denser the column of relics the 
super-high-energy neutrinos must traverse. In particular, 
the appearance of dips at energies much lower that we 
expect points to early — presumably nonacceleration — 
sources, that could give us insight into fundamental 
physics at early times and high energy scales. On the 
other hand, integration over a longer range in redshift 
means more smearing and distortion of the absorption 
lines. 

If the lightest neutrino mass is small, then Fermi 
smearing due to the thermal motion of the relics sets an 
effective lower bound on the neutrino mass, as reflected in 
the absorption lines. That moves the absorption lines to 
lower, more accessible, energies, but reduces the power of 
absorption-line spectroscopy to distinguish the neutrino 
mass eigenstates. Once the age of the cosmic-neutrino 
sources is clear, analyzing the thermal smearing may pro- 
vide at least a rough measurement of the current tem- 
perature of the relic neutrinos. In the best (imaginable) 
circumstance, discriminating the interactions of electron, 
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muon, and tau neutrinos might enable future experiments 
to determine, or verify, the flavor content of the mass 
eigenstates. 

It appears relatively secure to conclude that future 
observations of the v^jv^ ratio at the first absorption 



dip can decide whether the neutrino mass hierarchy 
is normal or inverted. Accelerator-based long-baseline 
neutrino-oscillation experiments should be the first to 
measure the sign of the "atmospheric" mass-squared dif- 
ference. This is not necessarily a losing proposition 
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for cosmic-neutrino absorption spectroscopy: the more 
information — neutrino masses, flavor composition, etc. — 
other experiments provide to neutrino observatories, the 
more perceptive cosmic-neutrino absorption spectroscopy 
might become about the thermal history of the universe. 

The experiments we describe in this study will not be 
done very soon, and their interpretation is likely to re- 
quire many waves of observation and analysis. Neverthe- 
less, they offer the possibility to establish the existence 
of another relic from the big bang and, conceivably, they 
may open a window on periods of the thermal history of 
the universe not readily accessible by other means. 
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